B2045+265: A New Four-Image Gravitational Lens from CLASSf] 
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ABSTRACT 

We have discovered a new gravitational lens in the Cosmic Lens All-Sky Survey 
(CLASS). The lens B2045+265 is a four-image system with a maximum separation 
of 179. A fifth radio component is detected, but its radio spectrum and its positional 
coincidence with infrared emission from the lensing galaxy strongly suggests that it is the 
radio core of the lensing galaxy. This implies that the B2045+265 lens system consists of 
a flat-spectrum radio source being lensed by another flat-spectrum radio source. Infrared 
images taken with the Hubble Space Telescope and the Keck I Telescope detect the lensed 
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and the National Aeronautics and Space Administration. The Keck Observatory was made possible by the generous financial support 
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images of the background source and the lensing galaxy. The lensed images have relative 
positions and flux densities that are consistent with those seen at radio wavelengths. The 
lensing galaxy has magnitudes of J = 19.2*™, mpimw = 18. 8 m and K = 17.6 m in a l'/9 
diameter aperture, which corresponds to the size of the Einstein ring of the lens. Spectra 
of the system taken with the Keck I Telescope reveal a lens redshift of Z£ = 0.8673 and 
a source redshift of z s = 1.28. The lens spectrum is typical of a Sa galaxy. The image 
splitting and system redshifts imply that the projected mass inside the Einstein radius 
of the lensing galaxy is M E = 4.7 x 10 11 /i _1 Mq . An estimate of the light emitted inside 
the Einstein radius from the K magnitude gives a mass-to-light ratio in the rest frame 
B band of (M/Lb)e = 20 h(M / Lb)q. Both the mass and mass-to-light ratio are higher 
than what is seen in nearby Sa galaxies. In fact, the implied rotation velocity for the 
lensing galaxy is two to three times higher than what is seen in nearby spirals. The large 
projected mass inside the Einstein ring radius may be the result of a significant amount 
of dark matter in the system, perhaps from a compact group of galaxies associated with 
the primary lensing galaxy; however, it may also arise from a misidentification of the 
source redshift. A simple model of the gravitational potential of the lens reproduces the 
image positions well, but further modeling is required to satisfy the constraints from the 
image flux density ratios. With further observations and modeling, this lens may yield 
an estimate of H . 

Subject headings: distance scale — galaxies: distances and redshifts — gravitational 
lensing — quasars: individual (B2045+265) 



2 



Introduction 



The Cos mic Lens All-Sky Survey (CLASS; |Mvers 



t aL 1998 ) is a large-scale survey for gravitational 
lenses among flat-spectrum radio sources. The pri- 
mary goals of CLASS are to find lenses which may 
be suitable for determinations of the Hubble Con- 



stant, Hp (e.g., Rcfsdal 1964; Blandford fc Narayan 



1 992 ), and to study the lensing rate in a large, ho- 
mogeneous sample in order to place limits on the cos- 
mological constan t , A (e.g., Turner, Ostrikcr fe Gott 
1984]; [Turner 1990|; |Fukugita, Futamase fe Kasai 199C| 



Fukugita fe Turner 199l| ). Over 12,000 sources have 
been observed with the VLA in three sessions - the 
first in the spring of 1994 (CLASS 1; -3300 sources), 
the second in the summer of 1995 (CLASS 2; -4500 
sources) and the third in the spring of 1998 (CLASS 3; 
—5000 sources). The vast majority of flat-spectrum 
radio sources are dominated by emission from a single 
compact core; all CLASS sources containing multiple 
compact components are selected as gravitational lens 
candidates, amounting to 50 - 100 candidates in each 
phase of the survey. These are followed up with high 
resolution radio imaging using MERLIN; the best sur- 
viving candidates are then imaged using the VLBA. 
During this procedure, the majority of VLA candi- 
dates are rejected based on surface brightness and 
morphology criteria. The candidates that survive the 
radio filters are then investigated further with optical 
and/or infrared imaging, spectroscopy and more ra- 
dio imaging. By these means we have discovered 11 
new lenses in the survey, in addition to the one re- 
ported here. We are investigating an additional —25 
promising lens candidates. 

We use Ho = lOO/i km s _1 Mpc" 1 and, except 
where noted, assume qo = 0.5 throughout this paper. 

2. VLA Observations 

The lens system B2045+265 (GB6 J2047+2643) 
was observed on 1995 September 02 as part of CLASS 2. 
The observations were made at 8.4 GHz with the VLA 
in A configuration, giving a resolution of — 0'.'25. Ob- 
servation and data reduction techniques for the survey 
will be discussed in Myers et al. (1998). The survey 
image of B2045+265 shows four components, with a 
possible detection of a fifth, weak component. The 



system has a standard lens geometry (e.g., Bland 
ford fc Narayan 1992 ) and is similar in appearance 



December 31. Observations were made in A configu- 
ration at 1.4, 4.9, 8.5 and 14.9 GHz. Details of the 
observations are given in Table [ij The source 3C 286 
was used as a flux calibrator. Phase calibrators were 
selected from the VLA calibrator list. The data were 
calibrated using standard AIPS routines, and maps 
were made using DIFMAP (Shepherd, Pearson & 



Taylor 1994). At the frequencies above 1.4 GHz, the 



to B1422+231 (Patnaik et al. 1992) 



The source was re-observed with the VLA on 1996 



source can be characterized as five point components, 
with no significant extended structure. Hence, the 
maps were made by fitting five point components to 
the data and then repeating a cycle of model-fitting 
and phase self-calibration. The self-calibration time 
scale was set to the length of the individual scans at 
each frequency. Both the positions and the flux den- 
sities were allowed to vary in the model fitting. At 
1.4 GHz the beam size is large compared to the com- 
ponent separation, so the model fitting was done by 
fixing the component positions at their 8.5 GHz val- 
ues and varying only the component flux densities. In 
addition, there were confusing sources in the 1.4 GHz 
map that had to be included as components in the 
model. 

The final maps are shown in Figures @ - @. In all 
cases, the morphology seen in the discovery map is du- 
plicated, with five distinct components clearly present 
in the high-dynamic-range 5 and 8.5 GHz maps. 
Component positions and flux densities were obtained 
using the model fitting procedures in DIFMAP. The 
brightest component (A) has an 8.5 GHz flux density 
of 16.55 mJy and is located at 20:47:20.29, +26:44:02.7 
(J2000); the positions of the other components rela- 
tive to component A are given in Table |^ (see Figure 
|l| for component labels). The positional uncertain- 
ties for the different components are estimated as the 
beam size at each frequency divided by the signal-to- 
noise ratio of the component flux densities; the posi- 
tional errors in the 8.5 GHz map are given in Table ^[ 
The largest separation in the system is l'/9. At all 
frequencies, five point components are adequate to fit 
the data. No sign of extended emission is seen, even 
in the deep 8.5 GHz map. The final component flux 
densities and RMS noise levels for the maps are given 
in Table [| In all cases the noise levels are within 5% 
of the thermal noise expected in the maps. 

The 8.5 GHz data were also analyzed to search 
for polarized emission. The polarization calibration 
was carried out in AIPS, using 3C 286 as the cali- 
brator. The phase calibrator 2115+295 was observed 
over a range of parallactic angles and was used to 
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Table 1 
Radio Observations 
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Fig. 1. — December 1996 8.5 GHz map. The contours 
are (-2.5, 2.5, 5, 10, 20, 40, 80, 160, 320, 640) times 
the RMS noise level of 0.0145 mJy/beam. Map made 
by fitting point source components, with flux densities 
listed in Table ||, to the (it, v) data, and restoring with 
a 0'.'27x0'.'23 restoring beam. 
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Fig. 2. — December 1996 4.9 GHz map. The contours 
are (-2.5, 2.5, 5, 10, 20, 40, 80, 160, 320) times the 
RMS noise level of 0.0621 mJy/beam. Map made by 
fitting point source components, with flux densities 
listed in Tabic 0, to the (u,v) data, and restoring 
with a 0'.'44x0'.'41 restoring beam. 
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Fig. 3. — December 1996 15 GHz map. The con- 
tours are (-2.5, 2.5, 5, 10, 20, 40) times the RMS 
noise level of 0.139 mJy/beam. Map made by fitting 
point source components, with flux densities listed 
in Table [|, to the (u, v) data, and restoring with a 
0'.'14x0'/13 restoring beam. 



Fig. 4. — December 1996 1.4 GHz map. The contours 
are (-2.5, 2.5, 5, 10, 20, 40, 80, 160, 320) times the 
RMS noise level of 0.0784 mJy/beam. Map made by 
fitting point source components, with flux densities 
listed in Table |^, to the (u, v) data, and restoring 
with a l'.'5xl'.'4 restoring beam. 



determine the instrumental polarization. The final 
maps are consistent with no polarized emission above 
56/iJy/beam (0.34% of the peak unpolarized inten- 
sity). 

Radio spectra for the five components are shown in 
Fig. ||. Components A - C have very similar spectra, 
with spectral indices of af® ~ —0.2 and a\% ~ —0.6 
(S v ocu a ). The spectra of components D and E differ 
from those of the three brighter components. In order 
to determine which, if either, of these two components 
is the counter image to components A, B, and C, flux 
density ratios are computed for each component with 
respect to component A. By computing flux density 
ratios at each frequency, it is possible to avoid uncer- 
tainties due to errors in the absolute flux calibration. 
The resulting curves (Fig. ^) , normalized to their val- 
ues at 8.5 GHz, clearly indicate that the spectrum 
of component E differs from those of the other four 
components. 

3. MERLIN and VLB A Observations 

A 5 GHz MERLIN observation of B2045+265 was 
made on 1997 November 30, with approximately 
12.5 hr total integration on source. The final map 
(Fig. |t]) has an RMS noise level of 0.146 mJy/beam 
and an angular resolution of 60 milliarcsec (mas). The 
data are well fitted by 5 point sources, with flux den- 
sities given in Table |3|. However, there are some in- 
dications of additional emission around components 
A and C in the residual map, at the 2 - 3cr level. 
The excess emission may indicate that the images 
are slightly resolved. If mas-scale structure can be 
detected in the images, it can be used to put con- 
straints on the lens model. The transformation ma- 
trix between corresponding positions in the resolved 
images provides crucial limits on the lensing poten- 
tial by fixing its second derivatives at those points 
(e.g., Gorcnstcin ct al. 1988| ). Observations with high 
angular resolution and dynamic range are needed to 
search for such structure. A 5 GHz VLBA snapshot 
(35 min total integration on source) of this system 
was made on 1995 November 12. In this observa- 
tion, which has an angular resolution of ~1 mas, the 
three brightest components are detected. The flux 
density of each component is only ^60% of its flux 
density in the 5 GHz VLA observations. Once again, 
this may indicate the presence of extended mas-scale 
structure in this source which could be detected in 
more sensitive observations. Deep VLBA observa- 




Frequency (GHz) 

Fig. 5. — Radio spectra of the five components in 
the B2045+265 system. The error bars are 3 a errors, 
based on the RMS noise in the maps. 




Frequency (GHz) 

Fig. 6. — Flux density ratios for components B (open 
squares), C (filled triangles), D (open triangles) and 
E (open circles) with respect to component A. All 
curves have been normalized to their 8.5 GHz point. 
Error bars are calculated by propagating the RMS 
noise in the maps at each frequency. 
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Table 2 
Component Positions 



Component Aa a A(5 a a x h 



A 


0.000 


0.000 




B 


-0.134 


-0.248 


0.001 


C 


-0.288 


-0.789 


0.001 


D 


+ 1.628 


-1.007 


0.006 


E 


+ 1.121 


-0.824 


0.005 



a Positions relative to component A are 
taken from the 1996 Dec 31 8.5 GHz map. 
Component A is at 20:47:20.29, +26:44:02.7 
(J2000) 

b Uncertainties in relative positions calcu- 
lated using the 1996 Dec 31 8.5 GHz data and 
assuming that the uncertainties in the compo- 
nent flux densities are three times the RMS 
noise level in the map. 
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Component Flux Densities 
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Table 4 

Component Flux Density Ratios 



v (GHz) Date B/A C/A D/A E/A 

1.4 GHz 1996 Dec 31 0.54 ± 0.009 0.69 ± 0.010 0.097 ± 0.008 0.075 ± 0.008 
4.9 GHz 1996 Dec 31 0.54 ± 0.009 0.65 ± 0.009 0.11 ± 0.008 0.13 ± 0.008 

8.5 GHz 1995 Sep 02 0.51 ± 0.051 0.80 ± 0.059 0.13 ± 0.046 0.099 ± 0.046 
8.5 GHz 1996 Dec 31 0.54 ± 0.003 0.66 ± 0.003 0.11 ± 0.003 0.14 ± 0.003 

14.9 GHz 1995 Sep 11 0.59 ± 0.072 0.76 ± 0.078 0.079 ± 0.062 0.15 ± 0.062 

14.9 GHz 1996 Dec 31 0.54 ± 0.044 0.67 ± 0.047 0.098 ± 0.046 0.18 ± 0.063 
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tions of B2045+265 have now been scheduled which 
may produce maps with sufficient dynamic range to 
detect extended structure in the images. 

4. Infrared Imaging 
4.1. NIRC 

The system was imaged on f996 July 3f using 



C\2 



the Near Infrared Camera (NIRC; Matthews & Soifer 
1994) on the W. M. Keck I Telescope. Images were 
taken in both J and K bands, with 45 one- minute 
exposures in K and 27 one-minute exposures in J. In 
order to estimate the point-spread function (PSF), 18 
exposures of a star at a distance of 1.47 arcmin were 
interleaved with the exposures on the lens in each 
band. The seeing was 0"45-0"90 in K and 0'.'75-l"15 
in J during the observations. 

The dark current level was subtracted from each 
image, and then sky-subtraction and gain correction 
were performed. The sky and gain frames used for 
each image were constructed from images observed 
directly before and after it. The individual frames 
were aligned by centroiding on a star which appeared 
in each frame. For the highest sensitivity in the final 
images, all 45 frames in K and all 27 frames in J 
were combined into mosaics. The K band mosaic is 
shown in Fig. || and close-ups of the lens system in 
the two bands are shown in Figs. || and [l(]. The two 
main features in the images of the lens system are a 
short arc, corresponding to the three brightest radio 
images, and the lensing galaxy, which is located l'/2 
from the arc. Spectroscopy has shown that the object 
seen 2" to the left of the lens system in Fig. || is a star. 

4.2. NICMOS 

B2045+265 was observed with the Hubble Space 
Telescope using the Near Infrared Camera/Multi- 
Object Spectrometer (NICMOS) on 1997 July 14. 
The NIC-1 camera was used, which has a pixel scale 
of 43 mas. Two exposures were taken giving a total 
exposure time of 2624 sec. The images were subjected 
to the standard NICMOS calibration pipeline involv- 
ing bias and dark current subtraction, linearity and 
flat-field correction, photometric calibration and cos- 
mic ray identification and removal. The final image 
(Fig. [ll]) clearly shows the three brightest lensed im- 
ages and the lensing galaxy. The relative positions of 
components A, B, and C match those seen at radio 
wavelengths to within O'.'Ol. In addition, the galaxy 
location matches the position of radio component E 
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Fig. 7.— MERLIN map. Components A, B, and C 
may be slightly resolved. 
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Fig. 8.— NIRC K band mosaic of the B2045+265 
field. The lens system is marked with the crosshairs, 
with the lensing galaxy to the left and an arc of emis- 
sion from the background source to the right. The 
object 2" to the left of the galaxy is a star. The cir- 
cled objects are extended and may be galaxies in a 
group at the redshift of the lens (see § 7.3). 
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Fig. 9.— NIRC K band image of B2045+265. The 
lensing galaxy is to the left, perhaps with some emis- 
sion being contributed from component D. An arc 
consisting of emission from components A, B and C 
is to the right. 



o 

0) 
O 



CY2 




Fig. 10. 



A x (arcsec) 
NIRC J band image of B2045+265. 



to within the errors. There appears to be low surface 
brightness emission between images B and C, indi- 
cating the possibility that extended optical emission 
from the background source is being lensed into an 
arc-like structure. The expected position of compo- 
nent D is marked in the figure. 

4.3. Photometry 

Because the NICMOS image has high spatial res- 
olution, it is possible to compute magnitudes for 
the individual images of the background source in 
the F160W bandpass (roughly corresponding to the 
ground-based H band). Magnitudes were calculated 
using both the DAOPHOT package ( |Stctson 1987|) 
in IRAFpl an d the SExtractor package (Bertin & 
Arnouts 1996). First, the magnitudes of the three 



bright lensed images were calculated using 0'.'26 diam- 
eter apertures. Next, an empirical PSF, constructed 
from radial profiles of stars in the field, was fitted 
to the emission from the lensed images and the star 
near the lensing galaxy. The scaled PSFs were then 
subtracted from the data and the lens magnitude was 
calculated using a l'/9 diameter aperture. This aper- 
ture was chosen to match the largest image separation 
in the system, which is approximately equal to twice 
the Einstein radius of the lens. The aperture diame- 
ter corresponds to 7.9 kpc at the redshift of the 
lens. The PHOTFNU header card was used to con- 
vert count rates in the image into flux densities into 
Janskys and then a Vega zero-point of 1087 Jy was as- 
sumed to get the F160W magnitudes. For comparison 
with the NIRC photometry (see below), the magni- 
tude of the "arc" defined by the three lensed images 
of the background source was also computed. The 
DAOPHOT and SExtractor packages produced mag- 
nitudes that were consistent within the errors of the 
sky determination. The final magnitudes are given 
in Table ^. It should be possible to search for ex- 
tinction caused by the lensing galaxy by comparing 
the F160W and radio flux density ratios of the lensed 
images. However, the uncertainties in the NICMOS 
component magnitudes are large (~ 0.2 m ) and thus 
the F160W flux density ratios, B/A = 0.69 and C/A 
= 0.76, are consistent with the radio flux density ra- 
tios at the 1 a level. The upper limit on the brightness 



5 IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, which are 
operated by the Association of Universities for Research in As- 
tronomy under cooperative agreement with the National Sci- 
ence Foundation. 
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of a point-source at the location of component D is 
mFwow > 22 m . This is not surprising since compo- 
nent D should be 2.5 m fainter than A. 

It is considerably more difficult to do photometry 
on the NIRC images because the seeing disk is large 
compared to the component separations. Instead of 
finding magnitudes of the individual lensed images, 
the magnitude of the "arc" was calculated by cen- 
tering a rectangular aperture of width 0'.'9 and height 
1'.'8 on the "arc." The sky was estimated from regions 
directly to the north, west and south of the aperture. 
Once again, the DAOPHOT package was used to sub- 
tract the emission from the lensed images and the star 
from the data before the lens galaxy magnitude was 
calculated. The photometric zero points were esti- 
mated by observing the i nfrared standard star SJ9184 
(K = 11.82, J = 12.18; |Pcrsson ct al. 1998| ), which 
was observed in each band prior to the lens and PSF 
star exposures. The zero points in the two observed 
bands were K = 22.44, and J = 22.92. The J and 
K magnitudes of the lensing galaxy and "arc" are 
given m Table |. 

5. Optical Spectroscopy 

Spectra of B2045+265 were taken with the Low 



Resolution Imaging Spectrograph (LRIS; |Oke et al 



1995[ ) on the Keck I Telescope on 1996 June 18-19. 
A 1" longslit and the 300 gr/mm grating were used, 
giving a pixel scale of 2.44 A/pix. The slit was po- 
sitioned at P.A. = 112°, in order to have the highest 
probability of spatially separating the emission of the 
background source from that of the lensing galaxy. 
Seven 1500 sec exposures were taken, four on the first 
night and three on the second night. 

The spectra were reduced using standard IRAF 
tasks. The average seeing over the two nights was 
sufficiently good to enable the extraction of separate 
spectra at the positions of the lensing galaxy and 
the brightest source emission. The spatial separa- 
tion between lens and source spectra is l'/3, matching 
the separation between the lens and the three bright 
lensed images seen in the NICMOS image. The fi- 
nal spectra have a wavelength range AA4309-936lA. 
Flux calibration was performed using the standard 
star Fcigc 110 ( pkc 1990[ ) which was observed at 
the end of each night. An atmospheric absorption 
template was created by fitting a power law to the 
flux-calibrated spectrum of the BL Lac object 2155- 
304. This template was used to remove atmospheric 



features from the B2045+265 spectra. The individ- 
ual exposures were weighted by the squares of their 
signal-to-noise ratios and combined to create the final 
lens and source spectra shown in Figs. 12] and 13. 



The lens system was re-observed with LRIS on 
1996 September 15 with a setup similar to that used 
in the earlier observations, but with wavelength cov- 
erage AA2344-7398A. Two 2000 sec exposures were 
taken and the data were reduced using a method sim- 
ilar to that described above. In this case, however, the 
chip response was not removed from the final spec- 
trum. The seeing was worse than in the earlier obser- 
vations; as a result, the source and lens spectra could 
not be extracted separately and both source and lens 
features are seen in the final spectrum (Fig. 



14) 



The reduced spectrum of the lensing object is typ- 
ical of a Sa galaxy (e.g., Kennicutt 1992). It shows 
[O II] A 3727A emission at 6962A, Ca II K A3934A ab- 
sorption at 7345A, Ca II H A3968A absorption (pos- 
sibly blended with He A3970A) at 7411 A, and H77 and 
m XX 3835, 3797A absorption at 7159A and 7092A, 
respectively. These spectral features establish the lens 
redshift as zi = 0.8673±0.0005, where the uncertainty 
in the redshift is determined from the RMS scatter in 
the redshifts derived from the [O II], Ca II K, Hi] 
and H9 lines. The source spectrum shows a broad 
emission line at A = 6396A. No other broad emission 
features are seen in the spectrum, which ranges from 
4300 to nearly 9400A. We identify the line as Mg II 
A2800A, which implies z s = 1.28, because it is the 
only broad emission line that is so isolated in typical 



quasar spectra (e.g., Boyle 1990: Francis et al. 1991). 
The second set of spectra were taken in order to search 
for C III] A1909 emission, expected to fall at ~ 4350A. 
No line was seen at this position, although the sen- 
sitivity in this part of the spectrum is low and the 
presence of a weak line is not ruled out. Other possi- 
ble identifications for the observed emission line (e.g., 
H/3, C III], C IV, Lya) would imply the presence of 
strong emission lines in parts of the spectrum in which 
the signal-to-noise ratio is high, none of which we see. 
We conclude that the line is correctly identified and 
the source lies at redshift z s = 1.28±0.01, where the 
uncertainty in the redshift is estimated from uncer- 
tainties in finding the line centroid in both the emit- 
ted and the observed frames. 
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Table 5 

Component Data (Infrared Observations) 



Component Ao a A<5 a J ^Fl60W b K 



A c 0.00 0.00 ■■■ 20.6 

B c -0.13 -0.23 ■■■ 21.0 

C c -0.28 -0.78 •■■ 20.9 

Lens Galaxy d +1.10 -0.80 19.2 18.6 17.6 

"Arc" ■ ■ ■ ■ ■ ■ 19.6 18.7 17.3 



a Rclativc positions in arcseconds calculated with respect to 
component A after a —23° rotation was performed to match the 
radio image orientation. 

b F160W magnitudes computed assuming that the Vega zero 
point is 1087 Jy. 

c Magnitudes computed in a 0''26 diameter aperture. 

d Magnitudes computed in a l''9 diameter aperture. 



6. Source Variability 



o 

CD 
w 

V 

S-H 




Fig. 11.— NICMOS F160W image of B2045+265. 
The extended object near the center of the image is 
the lensing galaxy. The crosshairs mark the expected 
position of component D. 



In order for a lens system to be used to measure 
Hq, the background source must be variable so that 
time delays can be measured. We use two methods 
to search for evidence of variability in the background 
source in B2045+265. The first is to examine the com- 
ponent flux density ratios at different epochs. The 
ratios are not affected by errors in the absolute flux 
calibration and so any detected changes in the flux 
density ratios reflect actual changes in the flux den- 
sity of the background source. A comparison of the 
8.5 GHz observations made in September 1995 with 
those made in December 1996 shows small changes in 
the ratios of the component flux densities, but not at 
a significant level. 

The second method is to monitor the total 22 GHz 
flux density of the B2045+265 lens system. Daily ob- 
servations of the system were made with the 40 m tele- 
scope of the Owens Valley Radio Observatory during 
the periods 1996 May 18 - September 03 and 1996 
November 04 - 22. The data were calibrated using 
the CMBPROG package ( |Lcitch 1998[ ). A 15% ± 3% 
change in total flux density was observed in November 
1996 over a period of less than a week. Because sev- 
eral other sources being monitored showed no change 
in flux density during the same period, the observed 
change in flux density may indicate variability in the 
background source. However, it is also possible that 
the detected variability is due to changes in compo- 
nent E, which has a radio spectrum that is flatter 
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Fig. 12. — LRIS spectrum of the lensing galaxy from 
data taken on 1996 June 18-19. 




Fig. 13. — Spectrum of background source taken 
on 1996 June 18-19. There is some contamination 
from light from the lensing galaxy, which produces 
the [O II] emission line. 




IOD 5000 6000 7000 

X (A) 

Fig. 14. — Wavelength calibrated spectrum taken 
on 1996 Sept. 15. Both source and lens features are 
seen in the spectrum. The position of the expected 
C III] emission is marked. The signal drops to zero 
blueward of 4000A. 



that that of the lensed images. To determine which 
of the components is varying requires high angular 
resolution monitoring of the system with the VLA. 

7. Discussion 

7.1. The Nature of Component E 

The detection of five unresolved radio components 
in the B2045+265 system raises the possibility that 
component E is the fifth lensed image of the back- 
ground source. We believe, however, that compo- 
nent E is instead associated with the lensing galaxy. 
Although standard lensing models can produce five 
images in configurations similar to that seen in th e 
B2045+265 system (e.g., |Blandford fc Narayan 19~9^) , 
the central fifth image tends to be highly demagni- 
fied with respect to the other images. This behavior 
is not seen in B2045+265; in fact, component E is 
brighter than component D at high frequencies. In 
addition, the radio spectrum of component E (Figs. 
H and ||) differs sufficiently from the spectra of the 
other four components to suggest that component E 
is not related to the other components. Finally, the 
radio emission of component E is spatially coincident 
with the infrared emission from the lensing galaxy 
(see §4.2). We therefore conclude that component E is 
indeed associated with the lensing galaxy. It thus ap- 
pears that we have discovered a radio galaxy lensing a 
radio-loud quasar. Of previously known lens systems, 
only 2016+112 (Lawrence et al. 1984) may have radio 
emission associated with a lensing galaxy. However, 
recent observations of 2016+112 suggest that the ra- 
dio emission of its component CI may be lensed emis- 



sion from a (separate) background source (Garrett et 
al. 1996| ). 

The flat radio spectrum of component E in B2045+265 
is typical of an active galactic nucleus. If component 
E is, indeed, associated with the lensing galaxy, then 
the radio position of E gives an accurate location for 
the nucleus of the lensing galaxy. The location of 
the center of the lensing galaxy relative to the lensed 
images is a key component in models of the lensing 
potential (see § |7.2| ). An a priori knowledge of the 
relative position of the center of the lensing galaxy 
provides strong constraints on lens models. 

7.2. Preliminary Lens Model 

We model the lens potential as a singular isother- 
mal sphere potential with "mixed" shear (Kochanck 
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1991), in which the scaled lensing potential can be 
expressed as: 



ip = br + "fbr cos 2 ( 



r)- 



Locations in the image plane relative to the lens cen- 
ter are given by x — (r, 0) — (x, y). The critical radius 
of the isothermal sphere is given by 6, and 7 and 6* 7 
give the magnitude of the shear and its position angle. 
In the modeling process we have used an astronomi- 
cal rather than mathematical coordinate system such 
that x = Aa and y = AS. The polar coordinates 
are defined in the usual way, with r = (x 2 + y 2 ) 1 / 2 
and 9 = arctan(y/a;). This results in a left-handed 
coordinate system. 

Although the image positions are defined relative 
to the position of component A, they are treated 
as absolute positions. Thus, both the lens position 
(xi,ye) and source position (x s ,y s ) are included as 
parameters in the model. We assume that compo- 
nent E marks the nucleus of the lensing galaxy and 
fix the lens position at that location, leaving five vary- 
ing parameters and eight observational constraints (x 
and y for each of the four images). 

A downhill simplex routine ( Press et al. 1992] ) 



is 



used to minimize the difference between the observed 
and model positions, expressed as: 



X 



E 
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where the subscripts m and o refer to the model and 
observed values, respectively. The positional errors 
are defined to be the 8.5 GHz beam size divided by the 
sig nal-to-noise ratio of the component in the 8 . 5 GHz 



map. — To avoid biases introduced by the choice of 



starting values, we repeat the process with a grid 
of choices of initial conditions containing 3" values, 
where n is the number of varying parameters in the 
model fitting. The best-fit parameters for the model 
are given in Table ||, and the model is shown graph- 
ically in Fig. [15]. The observed image positions are 
recovered by this simple model of the lensing poten- 
tial, with RMS image displacements of <20 mas. 

Radio images are unaffected by microlensing, so it 
is possible to use the flux density ratios of the images 
to constrain the model of the lensing potential fur- 
ther. The observed flux density ratios, however, are 
difficult to reproduce with a simple model such as the 
one presented above. Although this simple model is 
adequate to recover the four image positions, it does 



not reproduce the relative magnifications of images 
A, B, C. The reason for this is easy to see. If, as 
we have assumed, the potential is smooth, then we 
can perform a Taylor expansion of it about the lo- 
cation of image B. If we approximate the images as 
collinear, then the highest order essential terms in this 
expansion are quadratic in the coordinate perpendic- 
ular to the A - C line and quartic parallel to this line. 
Higher order terms are, by assumption, ignorable. It 
is then possible to show that in the limit of small im- 
age separations there are several scaling laws includ- 
ing Sb ~ Sa + Sc, i-e., the flux of the central image 
approac hes the sum of the fluxes of the other tw o 
images (|Kovner fc Paczynski 198^ ; |Blandford 199C|) . 
This is clearly violated in this source, which indicates 
strongly that the potential is not smooth on the scale 
of the A - C image separation. 

In order to improve upon our model, we have in- 
troduced higher order terms in the potential such 
as might be produced by a local mass perturbation. 
Specifically we expand the potential around images A 
and C to third order along the A - C line and to fourth 
order around image B. We then impose continuity of 
the potential and its first two derivatives midway be- 
tween A and B and midway between B and C. There 
are then three observable quantities, the ratio of the 
image separations and the ratios of the image fluxes, 
that can be used to solve for the coefficients in these 
Taylor expansions. Performing this exercise leads to 
the conclusion that adding a small, positive mass per- 
turbation centered on a point located O'.'ll from B 
towards C suffices to recover the flux ratios as well 
as the image locations. (This perturbation might be 



caused by a spiral arm, for example (cf Mao & Schnei- 



der 1997)). The ratio of the (almost certainly unmea- 



surable) time delays Atg^/ ^bc is now found to be 
0.2; however, the delay AtsD should be quite robust 
to this perturbation within a given global model. In 
conclusion, none of the existing observations challenge 
the gravitational lens interpretation that we present 
here. 

7.3. Properties of the Lensing Galaxy 

Properties of the lensing galaxy can be derived by 
using the fact that it is acting as a gravitational lens. 
For example, the mass within the Einstein ring of a 
lens is (e.g., Blandford fc Narayan 1992 ): 



M E « 1.24X10 1 



D,D S 



13 



Table 6 
Lens Model 



Parameter 


Value 


b 


1"076 


7 


0.104 




69? 4 




+ I"l2 b 




-0''82 b 


X a 


+Q'.'72 




-0"67 



a The positions rela- 
tive to component A 
{x, y) and 8-y are defined 
in a left-handed C i Qor | d i- 
nate system. See § |7.2| 

b Hcld fixed during 
model fitting. 



C\2 



O 

w 
O 




Fig. 15. — Lens model showing critical curve (broken 
line), caustic curve (solid line), lens center (©), source 
position (filled diamond) and model image positions 
(*). The four circles are centered on the observed 
image positions. 



where 9e is the angular radius of the Einstein ring and 
Dg, D s , and Dg s are the angular diameter distances 
between observer and lens, observer and source and 
lens and source, respectively. The angular diameter 
distance between two objects at Zj and Z j (zi < Zj) is 
given by ( Carroll, Press, fc Turner 1992 ): 



A, 

where 



JT (l+*j)l n kl 1/2 ' 



\n k \^ 2 i(z uZ] )) , 



dz 



i [(i + z) 2 {i + n M z)-z{2 + z )n A }- 1 / 2 



and sinn(x) is defined as sinh(x) for flk < 0, x for 
Qfc = 0, or sin(ar) for f2& > 0, respectively. The Vt 
terms are defined in the standard fashion: 



n 



M 



8rrG 



A 
3Hf 



fit 



kc 2 



H Q H Q 



The integral X has an analytic solution for fl A = 0, 
and if Qui = 1 the angular diameter distance acquires 
the simple form 



2c 



1 



1 



1 



Dij H (i + Zj ) \(i + Zi y/* {i + z y/ 2 



For B2045+265, the lens redshift is secure. The 
source redshift, based on one emission line, is less cer- 
tain. For this reason, we express the physical quanti- 
ties derived in this section as functions of the source 
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redshift. We will also assume, for ease of discussion, 
an Einstein-de Sitter cosmology. The effects of vary- 
ing the source redshift and the cosmological model are 
presented in Table [?]. The angular diameter distances 
for this system are D? = 860/i~ 1 Mpc and 

D s _ 1.8673[(l + z s ) 1/2 - 1] 
~Dts ~ {l + z s y/ 2 - 1.8673 



Th e mass enclosed within the Einstein ring radius is 

calculated by assuming that the maximum image sep- 
aration is a measure of 20e', i.e., that 6*£=0'.'95. This 
assumption yields M E « 9.6 x 10 10 (D s /D is )h~ 1 M e 
for the B2045+265 system. Note that Me represents 
the total mass contained in the cylinder with pro- 
jected radius Re- The speed of a particle moving 
in a circular orbit (u c irc) at the Einstein ring radius 
is derived from Me by assuming that the mass distri- 
bution of the lensing galaxy is an isothermal sphere. 
A simple calculation is used to convert Me to the 
mass inside the sphere of radius Re- The value for 
follows directly from the mass inside the sphere. 



Th e luminosity of the lone ineido Re can bo oeti 

mated by converting its K magnitude to rest-frame 
Mb and My magnitudes. The fc-correction for the 
lens redshift and rest-frame (V — K) and (B — V) col- 
ors for typical Sa galaxies are taken from Poggianti 
(1997). The resulting luminosities are (Lv)e = 2.92x 
10 10 h- 2 L Q v and (L b )e = 2.36 x 10 w h- 2 L Q B . 
The rest-frame B band mass-to-light ratio becomes 
{M/L b )e = A.l(D s /D es )h(M/L B ) Q . 

Inserting the source redshift derived from the op- 
tical spectra (§||) yields a lens mass of Me ~ 4.7 x 
10 n h~ 1 M(? 3 within the Einstein ring radius (~ 4/i _1 kpc 
at the redshift of the lens). This mass implies a cir- 
cular velocity at R E of 570 km/sec, two to t hree 
ti mes the veloc ities seen in both nearby (e.g., Ru- 
bin et al. 1985 ) and more distant (0.1 < z < 1; 



Vogt et al. 1996 ) spiral galaxies. The Einstein ring 
mass-to-light ratios are (M/Lv)e = 16h(M/ L\/)q 
and (M/L b )e = 2O/i(M/i B ) . These mass-to-light 
ratios arc not unusual for lens sys tems (e.g., Fig. 7 of 
Kceton, Kochanek, fc Falco 1998|) , although they are 
higher than expected given the redshift of the lens 
and assumptions about the luminosity evolution of 
lens galaxies (Keeton et al. 1998). The large lensing 
mass required to produce the observed image separa- 
tion is a result of the relative closeness of the back- 
ground source to the lens in redshift space. There are 
several possible explanations for the unusually high 
projected mass. 



Compact groups of galaxies associated with the 
primary lensing galaxies have been discovered in re- 
cent observations of two lens systems (1115+080 and 
B1422+231; Kundic et al. 1997a,b; [Tonry 1998; ). We 
are conducting observations of the B2045+265 field 
to search for evidence of a group associated with the 
lensing galaxy (see also Fig. ^). The high contamina- 
tion rate by field galaxies (~85% at z 



1; 



Postman 



Lubin, & Oke 1998), however, makes high-redshift 



groups difficult to detect. If such a group exists and is 
gravitationally bound, its gravitational potential can 
enhance the image splitting by the lensing galaxy. In 
such a scenario, we would overestimate the mass of 
the lensing galaxy. 

Another explanation for the high projected mass 
inside the Einstein ring is that the lens may be one of 
a close pair of possibly merging galaxies. This situa- 
tion has been observed in the case of the CLASS lens 
1608+656. The initial observations revealed an Ein- 
stein ring projected mass of Me ~ 3.1 x lQ 11 h~ 1 M&, 
corre s ponding to v C j rc = 44 km/sec ( Myers et al. 
1995; Fassnacht et al. 1996). Later HST imaging 



showed that the lens is composed of a very close pair 
of galaxies, explaining the high mass and ellipticity re- 
quired for the 1608+656 system ( Jackson et al. 1997 ). 
We are unable to determine unequivocally whether 
two galaxies are present in the B2045+265 system 
because the sensitivity of our NICMOS observation 
is too low. More sensitive NICMOS and WFPC2 ob- 
servations of B2045+265 are needed to test the two- 
galaxy hypothesis. 

A third possibility is that the source redshift or 
the assumed cosmology (or both) are in error. We 
examine the results of varying the source redshift by 
assigning different identifications to the emission line 
seen in Fig. |l3|. We then re-calculate the physical 
parameters associated with the lens in three differ- 
ent cosmological models: (J1m,^a) = (1, 0), (0.2, 0), 
and (0.2, 0.8). We find that the effects of changing z s 
are larger than those produced by changing the cos- 
mological model. Of the alternative identifications of 
the emission line in the source spectrum, the choice 
of Lya produces lens galaxy properties that are the 
closest to those associated with normal nearby spirals 
(Table |). However, no sign of a Lya forest is seen on 
the short-wavelength side of the emission line, even 
after subtracting the estimated contribution of the 
lens galaxy from the source spectrum. The absence 
of other emission lines in the spectrum (§^) and the 
lack of a Lya forest argue against the identification of 
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the emission line as Lya. More sensitive spectroscopy 
of this system must be conducted to determine the 
source redshift unambiguously. 

7.4. Prospects for Measuring Hq 

A lens system can be used to estimate Hq by com- 
paring the time delays predicted by the lens model 
to the observed time delays. The time delays are a 
function of the lensing potential, the source and im- 
age positions, the redshifts of the source and lens, 
and the world model (through the angular diameter 
distances) : 



At, = (1 + zi) 



DiDs 
cD 



Is 



\Si - /3| 2 - ip{xi) 



th 



where Xi — (xi,yi) is the angular position of the i 
lensed image and (3 = (x s ,y s ) is the position of the 
source. The observable quantity is the difference in 
the time delays, Atjj = At, — Atj . The predicted de- 
lays are proportional to h^ 1 through the ratio of an- 
gular diameter distances, so if it is possible to measure 
delays, the ratio of the predicted and observed delays 
will give h. For the model given in § 7.2, the compo- 
nents are expected to vary in the order: C, A, B, D. 
The predicted delays between components A, B and C 
are small: AtAc — 6.5/i _1 hr and Atsc = 7.7/i _1 hr. 
The delays between each of the three bright com- 
ponents and component D are all on the order of 
142 h^ 1 d. It will be challenging to measure the short 
time delays since radio loud quasars typically do not 
vary significantly on those time scales. However, with 
high-sensitivity observations, it will be quite possible 
to measure the long delays in this system, if the back- 
ground source is variable. 

8. Summary and Future Work 

We have discovered a new four-image gravitational 
lens in the second phase of CLASS. The B2045+265 
lens system is unusual in several respects. First, it is, 
with the possible exception of the puzzling 2016+112 
system, the first known gravitational lens system in 
which both the background source and the lensing 
galaxy are radio sources. The radio emission from the 
lens, presumably from an active nucleus, can be used 
to locate the center of the lensing galaxy with high 
precision, which provides a strong constraint for lens 
models. Secondly, the flux density ratios of the three 
brightest lensed images are inconsistent with the pre- 
dictions of simple models of lensing potentials. This 



result may indicate that there is significant structure 
in the mass distribution. Finally, the projected mass 
(and mass-to-light ratio) inside the Einstein ring of 
the lens is unusually high for a Sa galaxy. This may 
imply that there is a significant amount of dark mat- 
ter associated with this galaxy, that the lens actually 
consists of a close pair of galaxies, or that the overall 
mass distribution of an associated group of galaxies 
is contributing to the lensing potential. Alternatively, 
the source redshift may be misidentified. 

There are intriguing hints that the images of the 
background source may be resolved by high angular 
resolution observations. We have obtained high dy- 
namic range VLBA observations which can be used 
to search for this possible extended mas-scale struc- 
ture. If such structure is detected, the transforma- 
tion matrices between the lensed images of the back- 
ground source provide strong constraints on the lens- 
ing model. These data may explain the unusual flux 
density ratios among the three bright images. The 
VLBA observations also will provide a more accu- 
rate position for the center of the lensing galaxy. In 
addition, more sensitive Hubble Space Telescope and 
spectroscopic observations can provide further de- 
tails about the lensing galaxy and background source. 
The goal of these observations is to develop a well- 
constrained model of the lensing potential. A pro- 
gram of monitoring can be used to search for time 
delays in the system. Any measured delays can be 
combined with the lens models to yield a measure- 
ment of Hq. 
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Table 7 
Lens Parameters 











D t 






M E 


(M/L B ) E 




Line ID 


Zs 


Qm 


Oa 


(h,- 1 Mpc) 


(h,- 1 Mpc) 


(ft -1 Mpc) 


(lO 11 fe _1 M ) 


(h(M/L B ) & ) 


(km/scc) 


Mg II 


1.28 


1.0 


0.0 


860 


890 


180 


4.7 


20 


570 






0.2 


0.0 


1000 


1100 


240 


5.3 


10 


560 






0.2 


0.8 


1200 


1300 


340 


5.1 


11 


500 


ClII] 


2.35 


1.0 


0.0 


860 


810 


330 


2.3 


10 


400 






0.2 


0.0 


1000 


1200 


470 


2.8 


9 


400 






0.2 


0.8 


1200 


1300 


650 


2.6 


6 


360 


CIV 


3.13 


1.0 


0.0 


860 


740 


350 


2.0 


8 


370 






0.2 


0.0 


1000 


1100 


530 


2.5 


8 


380 






0.2 


0.8 


1200 


1200 


700 


2.3 


5 


340 


Ly a 


4.26 


1.0 


0.0 


860 


640 


337 


1.8 


7 


360 






0.2 


0.0 


1000 


1100 


540 


2.3 


7 


370 






0.2 


0.8 


1200 


1100 


600 


2.1 


5 


330 
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